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ABSTRACT 

The first orbital solution for the spectroscopic pair in the multiple star system a Scor- 
pii, determined from measurements with the Sydney University Stellar Interferome- 
ter (SUSI), is presented. The primary component is of (3 Cephei variable type and 
has been one of the most intensively studied examples of its class. The orbital so- 
lution, when combined with radial velocity results found in the literature, yields a 
distance of 174l^| pc, which is consistent with, but more accurate than the Hipparcos 
value. For the primary component we determine 18.4±5.4Mq, —4.12 ±0.34 mag and 
12.7 ± 1.8 Rq for the mass, absolute visual magnitude and radius respectively. A Bl 
dwarf spectral type and luminosity class for the secondary is proposed from the mass 
determination of 11.9 ± 3.1 Mq and the estimated system age of lOMyr. 
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1 INTRODUCTION 

The class of variable stars with /3 Cephei as the prototype 
consists of massive, nonsupergiant stars whose low-order 
pressure and gravity mode p ulsations result in light, ra dial 
and/or line profile variations (jStankov fc Handleill2005l ). In 
the most recent catalogue of /3 Cephei type variables, the 
93 confirmed members have periods of strongest pulsation 
(some members display multi-mode pulsation behaviour) 
ranging from 1-8 hours and are of spectral type B0-B3 
(|Stankov fc Handlei|[2005h . 

Approximately 14% of the catalogued /3 Cephei stars 
are located in multiple star systems and therefore, with suf- 
ficient observational data, component masses can be deter- 
mined for comparison with those estimated from theory. 
Indeed this has been achie ved in the c ases of a Virginis 
jHerbison-Evans et allligTll ). (5 Cent a,uri llDavis et al.ll2005l : 
lAusseloos et al.ll2006l ) and A Scorpu l|Tango et al.ll2006D bv 
combining the results of spectroscopic and interferometric 
analysis. Recently, analysis of the eclipsing binary HD 92024 
has als o yielded the ma ss determination of the /3 Cephei pri- 
mary (|Frevhammer et al, 2005.). This last r esult is of par- 
ticular interest as Frevhammer et al] (|2005| ) note that the 
HD 92024 primary component spectrum closely resembles 
that of a Sco. 

As a 'classical' member (|Lesh fc Aizenmanlll978l ) of the 
/3 Cephei variable type stars, a Sco (HR 6084, HD 147165) 
has been one of the most intensively studied examples of 
this class and has also been used as a photometric standard 
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l|Viiapurkar fc Drilling! [l99^ ). lLutz fc Lutj ([1973) classified 
the system as a binary with B2 IV -I- B9.5 V components 
separated by 20" on the sky and a visual magnitude differ- 
ence of 5.31 mag. Speckle interferometry, lunar occultations 
and spectroscopy have since shown that the B2 IV com- 
ponent is in fact three stars - a double-lined spectroscopic 
pair and a 2.2 ma g fainter B7 tertiary, 0.4" distant from the 
spectroscopic pair (Beavers fc Cooklll980l : lEvans et al ] |l986l : 
iMathias et al.lll99ll ). The primary component has been iden- 
tified as a /3 Cephei type pulsator and has a classification of 
Bl HI. Therefore the system, as it stands in the literature, is 
quadruple: a spectroscopic pair, a tertiary component, and 
the fourth star is the fainter compon ent in the visu al com- 
mon proper motion pair ADS 10009 (Pisulski'1992'). 

Th e atmosphere of the P Cephei c omponent was anal- 
ysed bv lVander Linden fc Butled (|l988l ). who gave an effec- 
tive temperature of 26 150 ± 1 070 K with a pulsation cycle 
variation of 4 000 ± 2 000 K. These authors also gave an ex- 
cellent synopsis of early radial velocity and li n e profi le stud- 
ies of the spectroscopic pair. iMathias et al.l (|l99lh used a 
double-shock wave propagating in the stellar atmosphere to 
explain the observed stillstand in the radial velocity curves of 
the 656.3 nm Ha, 658.3 nm CII and 455.3 nm Si HI lines (a 
period of about an hour wherein the radial velocity remains 
relatively constant). The double-shock wave model was also 
successful in describing similar chara cteristics of two othe r 
(5 Cephei sta rs: BW Vul an d 12 L ac (|Mathias et al.lll992l '). 
The study of IMathias et al.l (|l99ll) was also the first to de- 
tect the spectral lines of the companion in the spectroscopic 
pair. However, neither the spectral type, luminosity class 
nor an esti mate of the fiux ratio (secondary/primary) was 
presented in IMathias et al](|l99ll ). a Sco also exhibits a con- 
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Table 1. Adopted parameters of reference stars used during ob- 
servations. The angular diameters (and associated error) were 
estimated from an intrinsic colour interpolation (and spread in 
data) of measu rements made with the Narr abri Stellar Intensity 
Interferometer l lHanburv Brown et al.|[l974l ). 



HR 


Name 


Spectral 


V 


UD Diameter 


Separation 






Type 




(mas) 


from cr Sco 


5928 


p Sco 


B2IV/V 


3.88 


0.25 ±0.03 


6? 50 


5993 


Li)^ Sco 


BlVp 


3.96 


0.26 ±0.04 


5? 93 


6153 


u) Oph 


A7Vp 


4.45 


0.50 ±0.02 


4? 83 


6165 


T Sco 


BOV 


2.82 


0.31 ±0.06 


4? 20 



spicuous Van Hoof eflFect - a small phase lag of the radial 
veloci ty of the hydrogen lines relative to that of all oth er 
lines (IVander Linden fc Butle^l 19881 : iMathias et al■lll991^ ■ 

The most re c ent a n alyses of the spectroscopic pair 
bv iMathias et al.l (|l99ll ). IChapellier fc Valtieij (|l992l ) and 
iPigulskil (|l99i r have produced or bital per i ods o f 33.012, 
33.011 and 33.012 days respectively. IPigulskil (Il992l ) has also 
concluded that the variations in the main pulsation period 
(increase in the first half of the 1900s, decrease after 1960, 
and again an increase from about 1984) are due to a combi- 
nation of evolutionary and light-time efi'ects. 

Hereafter, the following naming convention will be fol- 
lowed. The primary and secondary refer to the spectroscopic 
pair. The tertiary refers to the star that is approximately 
2.2 mag fainter and separated by 0.4" from the spectro- 
scopic pair. The fourth and final star will be referred to 
as the distant companion. The term o Sco will refer to the 
spectroscopic pair unless explicitly stated otherwise. 

In Section [3] we provide the first complete orbital so- 
lution for G Sco determined by long-baseline optical inter- 
ferometry. In combination with the only double-lined radial 
velocity measurements, the distance, spatial scale, mass and 
age of the components are quantified and compared to pre- 
vious estimates in Section |4] The details of our observations 
and a description of the parameter fitting procedure are de- 
scribed in Sections [5] and 13] respectively. 



2 OBSERVATIONS AND DATA REDUCTION 

Measurements of the squared visibility (i.e. the squared 
modulus of the normalised complex visibility) or were 
completed on a total of 31 nights using the Sydney Univer- 
sity Stellar Interferometer (SUSI. [Pavis et all 1 19991 ). Inter- 
ference fringes were recorded with the red beam-combining 
system using a filter with centre wavelength and full-width 
half-maximum of 700 nm and 80 nm resp ectively. This sys- 
tem was outlined by iTuthill et aU (jioM ) and is to be de- 
scribed in greater detail by Davis et al. (in preparation). 

The observatio ns and data reducti on followed the pro- 
cedure outlined by iNorth et al. I l|2007l ) using the adopted 
stellar parameters of the calibrator stars given in Table [1] 
A total of 262 estimations of were available for analysis 
(not including two observations that were removed due to 
possible tertiary contamination. Section [3. ip as summarised 
in Tabled 



3 ORBITAL SOLUTION 

The theoretical response of a two aperture interferome- 
ter to the combined light of a binary star is given by 
iHanburv Brown et al.lll970l ) 



2_V?+ l3^Vi + 2l3\Vi\\V2\ cos(27rb ■ p/\) 



(l + /3)2 ' 

where /3 < 1 is the brightness ratio of the two stars in the ob- 
served bandpass and Vi , V2 are the visibilities of the primary 
and secondary respectively. In the simplest case, stars can be 
modeled by a disc of uniform brightness with angular diam- 
eter Q (see Section [4.31 for the effects of limb-darkenening) . 
The component visibilities in equation ((T)) are then given by 



V 



2Ji{-K\b\e/\) 



(2) 



■K\b\e/\ ' 

where Ji is a first order Bessel function. The angular sepa- 
ration vector of the secondary with respect to the primary is 
given by p (measured east from north), b is the interferom- 
eter baseline vector projected onto the plane of the sky and 
A is the centre observing wavelength. The observed will 
vary throughout the night due to the orbital motion of the 
binary and Earth rotation of b. The Keplerian orbit of a bi- 
nary star, i.e. p as a function of time, can be parameterized 
with seven elements: the period P, eccentricity e, the longi- 
tude of periastron uj, epoch of periastron To, semi-major axis 
a, the longitude of ascending node Q, and the inclination i. 
When using two-aperture optical interferometry, the phase 
of the complex visibility is lost and hence, Q. and uj have an 
ambiguity of 180°. This is a direct result of the fact that 
component identities cannot be determined. Measurements 
of radial velocity can be used to remove the ambiguity of u 
but that of Q. remains. 



3.1 Effects of the Tertiary and Distant 
Companion 

The presence of a third or fourth component can greatly 
complicate analysis of interferometric data. The distant 
component is too far away and too faint to affect the in- 
terferometric data and can be neglected. Therefore, an in- 
vestigation into the effects of the tertiary star was conducted 
to validate the data and aid the final analysis. 

Firstly, the mean location of the tertiary, relative to 
the spectroscopic pair, was estimated by fitting a simple 
smoothed curve trajectory to the published vector separa- 
tions. The literature vector separations are given in Table [3] 
and are shown with the fitted curve in Fig. [T] The estimated 
vector separation for the mean observing date of the SUSI 
observations (B2006.357) is p = 486 mas and 9 — 238° where 
p and 9 are the magnitude and angle of the separation vec- 
tor. The magnitude difference of the tertiary with respect to 
the combined spectroscopic pair irradiance is approximately 
2.23 ±0.18 mag (mean an d standard deviation of the values 
from'Hartkopf et al. 2006 adjusted to our passband). Given 
the estimated separation (which is within SUSI's field-of- 
view), the irradiance of the tertiary star will be detected by 
SUSI and must be included in the analysis. 

The offset of the tertiary fringe packet in delay space, 
D, was calculated using 



D = |bjpcos(r?- 6*) , 



(3) 
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Table 2. Summary of observational data. The night of the observation is given in columns 1 and 2 as 
a calendar date and a mean MJD. Column 3 is the mean orbital phase calculated from the values in 
Table |4] The baseline and the mean projected baseline (in units of metres) are given in columns 4 and 5 
respectively. Reference stars and the the number of squared visibility measures for a night are listed in 
the last two columns. 



Date MJD Phase Nominal Projected Reference # 

Baseline Baseline Stars 



2005 


May 20 


53510.64 


0.121 


80 


79.80 


r 


Sco, 


p Sco, Lu Oph 


4 


2005 


May 24 


53514.58 


0.240 


80 


79.78 


T 


Sco, 


p Sco, LU Oph 


10 


2005 


May 25 


53515.59 


0.271 


80 


79.74 


T 


Sco, 


p Sco 


4 


2005 


July 19 


53570.46 


0.933 


80 


79.81 


T 


Sco, 


p Sco, LU Oph 


10 


2005 


July 20 


53571.43 


0.962 


80 


79.78 


T 


Sco, 


p Sco, Lu^ Sco 


11 


2005 


August 08 


53590.46 


0.539 


80 


79.89 


T 


Sco, 


iji Sco 


11 


2005 


August 11 


53593.40 


0.638 


80 


79.80 


T 


Sco, 


Ljl Sco 


8 


2005 


August 13 


53595.43 


0.689 


80 


79.85 


T 


Sco, 


LU^ Sco 


10 


2005 


August 16 


53598.37 


0.778 


80 


79.73 


T 


Sco, 


LU^ Sco 


2 


2006 


June 15 


53901.55 


0.963 


80 


79.76 


T 


Sco, 


p Sco, Lu^ Sco 


3 


2006 


June 16 


53902.56 


0.993 


80 


79.82 


T 


Sco, 


p Sco 


9 


2006 


June 17 


53903.50 


0.022 


80 


79.83 


T 


Sco, 


p Sco 


12 








u.uoo 




1 y.oo 




Sco 






2006 


June 24 


53910.52 


0.234 


80 


79.80 


T 


Sco 




11 


2006 


June 25 


53911.48 


0.264 


80 


79.78 


T 


Sco 




8 


2006 


June 26 


53912.50 


0.293 


80 


79.83 


T 


Sco 




3 


2006 


June 27 


53913.53 


0.326 


80 


79.88 


T 


Sco 




8 


2006 


June 28 


53914.48 


0.354 


80 


79.83 


T 


Sco 




12 


2006 


July 18 


53934.40 


0.958 


80 


79.82 


T 


Sco 




2 


2006 


July 19 


53935.43 


0.989 


80 


79.81 


T 


Sco 




11 


2006 


July 25 


53941.41 


0.170 


80 


79.78 


T 


Sco 




8 


2006 


August 01 


53948.41 


0.382 


80 


79.80 


T 


Sco, 


p Sco 


12 


2006 


August 03 


53950.40 


0.442 


80 


79.77 


T 


Sco, 


p Sco, Lu^ Sco 


8 


2006 


August 04 


53951.37 


0.472 


40 


39.88 


T 


Sco 




4 


2006 


August 07 


53954.42 


0.564 


40 


39.91 


T 


Sco, 


p Sco, Lu-^ Sco 


13 


2006 


August 08 


53955.45 


0.596 


30 


29.94 


T 


Sco, 


p Sco, uj-^ Sco 


12 


2006 


August 09 


53956.42 


0.625 


30 


29.93 


T 


Sco, 


p Sco, uj-^ Sco 


12 


2006 


August 10 


53957.42 


0.655 


30 


29.94 


T 


Sco, 


p Sco, Lu^ Sco 


12 


2006 


August 22 


53969.40 


0.018 


55 


54.92 


T 


Sco 




6 


2006 


August 26 


53973.41 


0.140 


55 


54.95 


T 


Sco 




6 


2006 


August 27 


53974.38 


0.169 


30 


29.93 


T 


Sco, 


p Sco, uj-^ Sco 


5 



assuming that the spectroscopic pair's fringe packet is lo- 
cated at the phase centre of SUSI. |b| and rj are the pro- 
jected baseline length and position angle respectively. It was 
found that only two observations had an offset within the 
centre half of the 140 fim observing scan. For approximately 
72% of the observations, the offset was sufficiently large to 
place the tertiary fringe packet entirely outside the observing 
scan. However, the SUSI data reduction pipeline windows 
the recorded scan about the detected peak fringe location 
l|lrelandll2005l ). Therefore, the tertiary fringe packet will not 
be within the 'data window' and should not affect the cal- 
culation of if both 

(i) the offset is sufficiently large, and 

(ii) the fringe packet of the spectroscopic pair has a 
greater amplitude than that of the tertiary (i.e. avoid the 
system 'locking' onto the tertiary). 

The first condition can easily be satisfied by the removal of 
all observations with an offset less than 30 pm - the outer 
quarter of a scan (beyond the edge of the window function). 
Hence, the two data points that were found within the cen- 
tre half of the fringe scan were rejected from the analysis. 
The second condition may not be satisfied when the primary 



and secondary fringe packets destructively interfere; i.e. at 
a minimum in the modulation of there is the possibil- 
ity of mistakenly measuring the interference pattern of the 
tertiary. Assuming the primary and secondary fringe packet 
destructive interference is absolute and the tertiary is com- 
pletely unresolved, then the expected < 0.013 will be 
due to the tertiary alone. As this is below the detection 
limit of SUSI, we can neglect such cases. Hence all remaining 
observations of a Sco now satisfy both conditions. The SUSI 
data pipeline will nevertheless be affected by the incoherent 
fiux from the tertiary component. The required adjustment 
to the fitted squared-visibility model is given in the next 
section. 



3.2 Fitting Procedure and Uncertainty 
Estimation 

The validity of equation ([T} is (strictly) only for observations 
of binary stars made with very narrow bandwidths. For real 
detection systems wide bandwidth effects can reduce the ob- 
served and for the scanning detection system of SUSI, 
the equivalent to equation ([T| giving for a binary star fo r 
the case of a wide spectral bandwidth is (|North et al.ll200^ ): 
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Table 3. Literature vector separations of the tertiary measured 
with respect to the spectroscopic pair. Methods are O: Occulta- 
tion; S: Speckle; H: Hipparcos; C: Coronograph 



Besselian 
Year 



Method: Reference 



(mas) (dcg) 



1972.556 

1976.471 

1977.4868 

1980.4792 

1980.4819 

1981.4567 

1981.4704 

1981.4730 

1983.4254 

1984.3783 

1986.243 

1987.2726 

1989.2275 

1989.3038 

1991.3140 

1993.3424 

1994.3509 

1997.2252 

1997.5174 

1997.5174 

1997.6157 

2000.3996 



322 
326 
353 
367 
367 
372 
369 
369 
377 
384 
392 
407 
416 
414 
428 
430 
441 
443 
465 
456 
447 
469 



296.8 
291.8 
285.4 
277.4 
277.9 
275.2 
277.1 
275.6 
272.6 
272.3 
268.5 
264.5 
261.4 
261.4 
258 
256 
253 
249.3 
242.9 
249.5 
248.8 
244.4 



O: Evans et al. (1986) 
S: Morgan et al. (1978) 
S: McAlister (1979) 
S: McA lister et al. (198 
S: McAUster et al. (198 
S: McAlister ct al. (198 
S: McAUster ct al. (198 



McAlister et al. 
McAlister et al. 
McAlister et al. 



Evans et al. (1986) 
McAlister et al. (1989) 
McAlister et al. (1990) 
McAlister et al. (1990) 
ESA (1997) 
Miura ct aL (1995 ) 
Mason (1996) 
Hor ch et al. (1999) 
Horch et al. (1999) 
Horch et al. (1999) 
Horch et al. (1999) 
C: Sha tskv fc Tokovinin (2002) 
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Figure 1. Literature vector separations of the tertiary measured 
with respect to the spectroscopic pair. Lunar occultations are 
given as triangles, speckle interferometry as diamonds, Hipparcos 
results as a square and the coronographic measurement as a cross. 
The plus symbol represents the location of the spectroscopic pair. 
The dotted line is a simple smooth trajectory fit to the data and 
the asterisk is the estimated mean position of the tertiary during 
SUSI observations. 



+ l3^Vi + 2(5r{4,)\V^\\V2\ cos(i/;) 



where 



r(%l)) = exp 



X2 
-^0 



32 In 2 



(4) 



(5) 



The spectral response is approximated as a Gaussian of cen- 
tre wavelength Ao with full- width half-maximum A A. The 
term r(-i/)) corresponds to the autocorrelation of the Gaus- 
sian envelope of the interference pattern and t/j = 2-Kb- p/Xo 
is defined for convenience. 

The measures of are contaminated by the irra- 
diance of the tertiary star (Section I3.1|l such that equa- 
tion (|4]) is no longer applicable. The brightness ratio I3 of 
the tertiary/(primary -I- secondary) is approximately I3 ~ 
0.128 ± 0.023 (Section [3T]). Adjusting equation (g]) for the 
contamination of the tertiary we obtain 

1^,2 ^ Vi^ + f^^Vj + 2/3| Vi \\V2\rW cos(V>) 

I [(l + l3)(l+^)]2 • W 

The term, I3, reduces the observed and hence the effect 
the tertiary component has on the calculated of the spec- 
troscopic pair is considered an extra incoherent source. The 
addition of this term will mainly affect the fitted component 
angular diameters and brightness ratio, leaving the orbital 
parameters relatively unaffected. As AA and I3 approach 
zero, i.e. narrow bandwidth observations of a simple binary 
star, then equation ((6)1 reduces to equation 

Initial values of the inclination and position angle of 
the ascending node were found by a coarse grid search of 
parameter space. The remaining orbital parameters were 
limited t o within three standard deviations of the values 
given by iMathias et aP (_1991). The inital angular diame- 
ter of the primary star was estimated from the Hipparcos 
distance and spectral type characteristics found in the liter- 
ature (Gox 2000). By physical arguments and inspection of 
the measured values, the secondary's angular diameter 
and brightness ratio were limited to the ranges 0.2-0.6 mas 
and 0.2-0.8 respectively. 

The final estimation of parameters was completed us- 
ing minimization as implemented by the Levenberg- 
Marquardt method to fit equation ((6| to the observed values 
of V^. When finding the minimum of the manifold, the 
inverse of the covariance matrix is calculated by the non- 
linear fitting program. The formal uncertainties of the fitted 
parameters are derived from the diagonal elements of this 
covariance matrix. As the visibility measurement errors may 
not strictly conform to a normal distribution and equation 
((6]) is non-linear, the formal uncer tainties may be un deresti- 
mates. Following the approach of iNorth et al.l (j2007l ). three 
uncertainty estimation methods were adopted to confirm 
the accuracy of the values derived from the covariance ma- 
trix. These methods: Monte Carlo, bootstrap and Markov 
chain Monte Carlo (MCMC) simulations are described by 
iNorth et all l|2007l ) and references therein. 

3.3 Results 

From preliminary analysis, a classification of Bl V was 
adopted for the secondary (see Section [4. 4p and its angular 
diameter was assumed to be 0.34 ± 0.04 mas (i.e. a radius of 

6.4 Rq - with a 12 per cent uncertainty - at the distance 
obtained from the dynamical parallax). The three param- 
eters 61, 62 and (3 are coupled and hence a change in one 
will affect the other two without significantly changing the 
orbital parameters. 

The best-fitting values of the model parameters are 
given in Table |4] and four nights data are shown in Fig. [2] 
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Table 4. Orbital and physical parameters of a Sco as found with SUSI and from the literature. The periods 
from the literature have been converted to MJD and/or adjusted to the same epoch. Note that the difference 
between the Heliocentric Julian Day (HJD) and a Julian Day is less than ten minutes which is much smaller 
than the uncertainties. 



Parameter 


Unit 


SUSI 


Mathias et al. ("1991) 


Pigulski (19921 


P 


days 


33.010 ± 0.002 


33.012 ± 0.002 


33.012 


e 




0.3220 ± 0.0012 


0.44 ±0.11 


0.40 ± 0.04 


To 


MJD 


34889.0 ± 1.0 


34888.9 ± 0.7' 


34888.0 ±0.4= 


UJ 


deg 


283 ± S'' 


299.1 ± 10.0 


287 ±6 


a" 


mas 


3.62 ±0.06 






Q 


deg 


104 ± 5 






i 


deg 


158.2 ±2.3 






Oi 


mas 


0.67 ±0.03 








mas 


0.34 ± 0.04" 






P 




0.48 ±0.02 







°- Adopted parameters; HJD-2400000.5; 498 X 33.012d added; <^ Ambiguity of 180°. 





Figure 2. Data from the nights of 2005 Aug 08 (top-left), 2006 Aug 08 (bottom- left), 2006 Jun 24 (top-right) 
and 2006 Jun 28 (bottom-right) where each data point represents a measure of with the associated formal 
error. The values of Table |4] have been been used to show the fitted model as a solid line. 
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-4 -2 2 4 

RA separation (mas) 

Figure 3. The relative orbit (solid line) of the secondary about 
the j3 Cephei star projected on the plane of the sky. The dotted 
line is the line-of-nodes, the open triangle signifies periastron and 
the grey squares are parts of the orbit that were observed with 
SUSI. The stellar components (filled circles) are drawn to scale 
using the values of Table |4] The secondary is shown at an orbital 
phase of 0.85. 



pearance. The Probability Density Functions produced by 
the MCMC simulation for the remaining (free) model pa- 
rameters were slightly non-Gaussian, most likely due to the 
intrinsic pulsations of the primary. Even though some model 
parameters produced (weakly) non-Gaussian distributions, 
the uncertainty values quoted in Table [3] are the standard 
deviations. As the MCMC simulation includes the uncer- 
tainties in the tertiary incoherent flux and angular diameter 
of the secondary, the values it produced form the basis of 
the final parameter uncertainties given in Table |4] because 
we believe they are the most realistic estimates of parameter 
uncertainties for our data set. 

The orbi t al pa ram eters found by the analyses of 
iMathias et"all (ligQll ) and lPigulskil (|l992l ), together with the 
final values determined from the SUSI data, are given in Ta- 
ble 21 The values for the period and time of periastron pas- 
sage are all in excellent agreement. There is some disagree- 
ment among the two remaining orbital parameters when 
considering the quoted uncertainties. However, all values are 
co nsistent at the t wo st andard deviation level. The analysis 
of iMathias et al.l l|l99ll ) is considered to be the most up- 
to-date due to the detection of the secondary spectral lines 
and the inclusion of the most recent (p ublished) data. There - 
fore, comparing only the SUSI and the IMathias et all (|l99ll ) 
results all parameters are consistent at the 1.1 standard de- 
viation level. 



with the predicted model overlaid as the solid curve. The 
projected orbit on the plane of the sky is shown in FigO 

The reduced of the fit was 3.58, implying that the 
measurement errors are underestimated by a factor of 1.89. 
Two possible effects were investigated. Firstly, the seeing 
conditions during some observations were poor and then 
only the brightest calibrator could be used, resulting in data 
of lower quality t han on other nights. While a non- linear 
seeing correction (|lrelanclll2006l ) was applied to all data as 
part of the data reduction, there could be some residual 
atmospheric effects. Secondly, the primary star, being of (3 
Cephei pulsator type, varies in diameter, temperature and 
apparent brightness. These properties will also affect the re- 
sults of fitting equation (|S} to the data and contribute to 
the somewhat large value of reduced and non-Gaussian 
parameter uncertainty distributions (discussed further be- 
low). However the quality and resolution of the SUSI data 
combined with the uncertainty in the literature values does 
not justify a more complicated model that includes the in- 
trinsic variability. Furthermore, the number and time base 
of measures is sufficient to average out any effects from 
the intrinsic variability of the primary. Hence, the angular 
diameter of the primary, along with the component bright- 
ness ratio, are treated as mean values by equation (|6]). 

The three different uncertainty estimation techniques 
(described in Section I3.2|l produced distributions of each 
free model parameter similar in appearance and approxi- 
mately centred on the best-fitting values. The Monte Carlo 
and bootstrap methods were set to each generate 10'^ syn- 
thetic data sets while the MCMC simulations completed 10^ 
iterations. The shape of the likelihood functions of the semi- 
major axis, eccentricity, inclination, uniform disc angular 
diameters and the brightness ratio were Gaussian in ap- 



4 SYSTEM AND PHYSICAL PARAMETERS 

IMathias et al.l l|l99ll) contains the only published semi- 
amplitudes of both the primary and secondary co mponents. 
The sm all discrepancy between the SUSI and Mathi as et al] 
(|l99ll ) eccentricity and longitude of periastron passage may 
affect the estimation of those physical parameters obtained 
by the combination of interferometric and spectroscopic re- 
sults. We note, however, that the inclination is close to 
180° and consequently, the uncertainty in sini will dominate 
the error budget (see for example Section [4. l|l . 



4.1 Distance 

The calculation of the dynamical parallax. 



requires the semi-major axis of the relative orbit in both lin- 
ear units (i.e. AU) and angular units. Using the interferomet- 
ric orbital values and the component semi-a mplitudes Ki = 
31.9 ± 1.3 kms"^ k K2 = 49.3 ± 7.3 kms~^ (|]Mathias et all 
ll99lD . the semi-major axes of the component orbits were 
calculated (in km) using 

0.1,2 = : — : . (8) 

TTsmi 

These values are given in Table [S] in units of AU. Combining 
these values of ai and 02 with the angular semi-major axis 
of the relative orbit in Table |4j the dynamical parallax was 
found to be tt^ = 5.76 ± 0.68 mas. The 12% uncertainty in 
the dynamical parallax is dominated by the contribution of 
the inclination - the uncertainty in sinz is 11%. 
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The distance to the system is therefore 174l^g pc, which 
is consistent within uncertainties with the Hippar cos value of 
225^5 pc. An analysis of the Hipparcos data bv lde Bruiind 
l|l999h produced a secular parallax of -Ksec = 5.27 ± 0.48 mas 
(or a distance of IQOljg pc) with which the dynamical par- 
allax (and consequent distance) is in good agreement. 

The me mbers of nearby OB as sociations have been de- 
termined by Ide Zeeuw et al.l l| 19991 ) from Hipparcos proper 
motion and parallax measurements. They concluded that 
a Sco is a secure member of the Upper Scorpius (US) sub- 
group within the Sco OB2 associatio n with a membership 
probability of 94%. Using fig. 5 of Ide Zeeuw et al.1 [l999l . 
the mean US parallax i s 6.9 mas with an est imated uncer- 
tainty of about 1.5 mas (|de Zeeuw et al.ll 19991 give a 'spread' 
value of 1.6 mas). Therefore the dynamical parallax posi- 
tions a Sco closer to the centre of US compared to Hipparcos 
and hence strengthens membership probability to the Upper 
Scorpius subgroup. 



4.2 Component Masses 

The mass of the primary and secondary can be extracted 
from the orbital solution by combining Kepler's third law, 



Ml +M2 



(ai + 02) 

p2 



(9) 



and the ratio of the component semi-major axes about the 
centre-of-gravity, 

(10) 

M2 ai 

The calculated masses are in solar units when the semi- 
major axes are given in astronomical units and the period 
in years. Using the values determined in Sections 13.31 and 
14. II the primary and secondary star star masses were found 
to be 18.4 ± 5.4 Mq and 11.9 ± 3.1 Mq respectively. Once 
again the uncertainty in the mass is dominated by the con- 
tribution of the inclination (propagated from the distance 
determination) . 

The expected r nass of a Bl III star i s 17.4 Mq (interpo- 
lated from values in ISchmidt-Kaleilll982l ). while the analysis 
of HD92024 ( a Bl III P Cephei star w ith a nearly identical 
spectrum) bv lFrevhammer et al.l (|2005l ) produced I5I4 -^0- 
Our determination of the primary star's mass in a Sco is in 
complete agreement with these results. The mass of the pri- 
mary is larger than the maj ority of members in the ca, talogue 
of Galactic f3 Cephei stars (jStankov fc Handleij|2005l) where 
a mass of 12 Mq is the norm. The uncertainty in the mass 
is large, precluding the conclusion that the (3 Cephei com- 
ponent in a Sco is one of the most massive examples of this 
type of pulsating star. 



4.3 (3 Cephei Component Radius 

The analysis of the interferometric data yielded the uniform 
disc angular diameter of the primary component. However, 
real stars are limb-darkened and corrections are required 
to find the 'true' angular diameter from the uniform disc 
value. Thes e corre ctions are discussed (and tabulated) in 
iDavis et al.l (|2000[ ) and are dependent on the star's effec- 
tive temperature, surface gravity, chemical composition and 



the wavelength at which the uniform disc diameter was de- 
termined. Assuming solar c hemical composition, for Teff = 
26 150 K and logg = 3.85 d Vander Linden fc Butleilll988h 
at an observing wavelengt h of 700 nm, the li mb-darkening 
correction factor given in iDavis et al.1 l|2000l) is 1.018 for 
the primary star. Using the dynamical parallax, the ra- 
dius can now be estimated to be 12.7 ± 1. 8 R|^, compat- 
ible w ith a Bl III star (interpolated) from ISchmidt-Kaleij 
lll982h and the analysis of the similar star HD92024 by 
iFrevhammer et al.1 (120051 ). 

4.4 Secondary Component 

Even though iMathias et al.l (|l99ll ) detected spectral fea- 
tures of the secondary star, no estimation of its spectral 
type or luminosity class appears in the literature. 

Initial analysis of the interferometric data included the 
secondary angular diameter as a fitting parameter and a 
uniform disc radius of 11 Rq was determined (but poorly 
co nstrained). Interpolat ing spectral type calibration tables 
in lSchmidt-Kaled (|l982l), stellar types Bl V & B3 III had a 
mass that was consistent with that of the secondary (Sec- 
tion |4]2]). However, the fitted uniform disc radius was too 
large for a Bl dwarf star (6.4 Rq) but was in accord for a 
B3 giant (10.8 Rq). On the other hand, single-star evolu- 
tionary tracks (Section 14. 6p predict vastly different ages for 
the two components and produce a more evolved state for 
a B3 giant secondary. Therefore a Bl dwarf classification 
was adopted and the interferometric data was reanalysed 
with the additional constraints on the likely radius of the 
secondary taken into consideration. The results with this 
new fitting produced plausible values (absolute magnitudes 
and primary component radius more in line with literature 
calibrations, and similar ages) and without any significant 
change in the orbital parameters or quality of the model-fit. 

Therefore, evolutionary tracks and the mass determined 
from the combination of interferometric and spectroscopic 
results imply a classification of Bl V for the secondary. An 
approximate effective temperature of T^ti ~ 25 400 ± 2 000 K 
(Lang 1991j) is henceforth adopted for the secondary com- 
ponent. 

4.5 Component Magnitudes 

The absolute visual magnitude of the system can be de- 
termined from the distance and the apparent visual mag- 
nitude. The mean and standard deviation of extinction 
values found in the literature are Ay = 1.23 ± 0.20 mag 
llClavton fc Hanson|[l99^ : IWegneij[200^ : ISartori et al1l2003l : 
Ide Bruiind I I999I ) and the the appar ent visual magnitud e 
of a Sco is = 2.88 ± 0.02 mag (| Johnson et all Il966l ). 
Combining these values with the distance from Section [4. II 
the absolute visual magnitude of the spectroscopic pair is 
Mv(Bl HI -FBI V) = -4.55 mag. 

The absolute visual magnitudes of the components 
can now be determined from the brightness ratio and 
Afy(Bl HI + Bl V). However, the estimation of the bright- 
ness ratio in Section 13.31 was made at a wavelength of 
700 nm and needs to be adjusted to the centre of the 
V band (550 nm). Using blackbody radiation curves for 
Tcff = 26 150 K and T^g = 25 400 K the adjustment fac- 
tor to the measured brightness ratio was found to be 0.996 
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Table 5. Physical parameters of tr Sco. The middle and bottom panels are the values for the primary 
and secondary respectively. 



Parameter 



Unit 



This Work 



Literature 



Reference 



ai 


AU 


0.25 ± 0.03 






0.2 


AU 


0.38 ± 0.07 






7r 


mas 


5.76 ± 0.68 


4.44 ±0.81 


ESA (1997) 




mas 




5.27 ±0.48 


de Bruiine fl999') 


distance 


pc 


174+23 


225+50 
190tl^ 


ESA (19971 
de Bruiine (1999) 


age 


Myr 


10 


5-14 
5 


Brown (1998') 
Preibisch et al. (20021 



M{B1 III) 


M© 


18.4 ±5.4 






Mv(Bl III) 


mag 


-4.12 ±0.34 






L{B1 III) 


L© 


(2.9 ±0.8) X 10* 






i?(Bl III) 


R© 


12.7 ± 1.8 






M{B1 V) 


Mq 


11.9 ±3.1 






Mv(Bl V) 


mag 


-3.32 ±0.34 






L(B1 V) 


L© 


(1.6 ±0.4) X 10" 







(i.e. a 0.004 mag change to the magnitude difference) and 
hence, the V band brightness ratio is 0.48 ± 0.02 (i.e. a 
0.8 mag visual magnitude diflerence). Hence, Mv(Bl III) = 
-4.12 ± 0.34 mag and Mv(Bl V) = -3.32 ± 0.34 mag (the 
uncertainty estimates include the uncertainties in the dis- 
tance, reddening and the component brightn ess ratio) . These 
valu es are in acc ord with the tabulations of IPanagial (Il973l ) 
and lLand l|l99ll ). 



4.6 Luminosity and Age 

Bolometric corr ections found in the literature f or a Bl gian t 
are -2.13 mag (|Panagialll973 ) and -2.43 mag |Lan3[l993). 
Adopting the mean value of —2.28 mag, the luminosity of 
the primary component is (2.9 ± 0.8) x IO^Lq. The sec- 
ondary star has a proposed classification of a Bl dwarf which 
has bolometric corre ctions of —2 .23 mag and —2.70 mag from 
iPanagial (|l973l ') and iLan^ (|l99lh respectively. The luminos- 
ity of the secondary is (1.6 ± 0.4) x 10* Lq using a mean 
bolometric correction of —2.47 mag. Both the primary and 
secondary luminosities are consi s tent ( withi n unc e rtaint ies) 
with the tabulations of lPanagial (|l973l l and lLand (119911 ). 

The ages of the component stars can no w be es t imate d 
using the single-star evolutionary models of IClare3 l|2004l ). 
The positions of the stars in the HR diagram (Fig. |4]) are 
marked with crosses, the lower-luminosity one being the sec- 
ondary star (with larger error bars refiecting the consider- 
able uncertainty in its efi'ective temperature). Isochrones of 
9, 10 and 11 Myrs have been calculated and shown as dotted 
lines with the dashed line corresponding to an age of 5 Myr. 
The evolutionary-model masses of both star!{3 are compati- 
ble with the determination in Section [4.21 even though the 
primary has a lower evolutionary-model mass than the mea- 
sured value. 



The initial erroneous classification of the secondary as a B3 
giant (Section l4.4l l produced a location in the bottom right-hand 
corner of Fig. |4] corresponding to an evolutionary-model mass of 
approximately 8Mq. 




4.55 4.50 4.45 4.40 4.55 4.30 4.25 
log Te„ (K) 

Figure 4. Single-star evolutionary tracks of IClare"3 l|2004h with 
calculated isochrones of 9, 10, 11 Myr given as dotted lines with 
the dashed line corresponding to the US group age of 5 Myr 
dPreibisch et al.l l2002h . The locations of the components are 
marked - the lower cross being the secondary. 

The age of the system is estimated to be 10 Myr and 
from the isochrones coeval formation of the two components 
is confirmed. The age range of the entire Sco 0B2 association 
(US + Uppe r Centaurus-L upus + Lower Centaurus-Crux) 
is 5-14 Myr (jBrownl Il998l ). consistent with the age deter- 
mined for a Sco. Ifowe ver, a recent exploration of the stellar 
population of US by Preibisch et al . (2002) suggests an age 
of 5 Myr for the group. They propose that the US members 
are the result of a shock wave passing through a molecular 
cloud 5-6 Myrs ago causing a burst in star formation and the 
shock wave being from a supernova in the Upper Centaurus- 
Lupus group about 12 Myr ago. Ifowe ver, a Sco is an ' o utlier ' 
in th e colour-magnitu de diagram of iPreibisch et"al] (I2OO2I ) 
(and lde Bru"iin3ll999l ) who explain the deviation to be a re- 
sult of its binary and pulsation characteristics. Hence there 
is an implication that a Sco may have formed prior to the 
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other members of the US group but the effect of binarity and 
observational uncertainties must be investigated further be- 
fore the evolutionary status of a Sco can be confirmed. 



5 SUMMARY 

The first complete orbital solution for cr Sco, based on inter- 
ferometric measurements with SUSI, is presented. In com- 
bination with the only double-lined radial velocity measure- 
ments, the dynamical parallax and distance to a Sco have 
been determined and shown to be consistent with previous 
estimates. Furthermore, the masses of the components have 
been determined which, in combination with evolutionary 
tracks, allows the first proposed classification of the sec- 
ondary as Bl V. 

Using calibrations found in the literature, the absolute 
visual magnitude, luminosities and mass of the component 
stars are found to be consistent with types Bl III and Bl V 
for the primary and secondary respectively. Furthermore, 
the radius of the primary, determined from the interfero- 
metric angular diameter and distance, is in accord with a 
Bl giant classification. 

The lOMyr age of the spectroscopic pair, e stima ted 
from the single-star evolutionary models of lClarej l|2004) . is 
in agreement with previous estimates of the Sco OB2 associ- 
ation as a whole (5-14 Myr) but not with other members of 
the Upper Scorpius group (5 Myr). Further information (in 
the form of improved component characteristics) is needed 
to either resolve this discrepancy or show that a Sco formed 
earlier than the remainder of the Upper Scorpius group. 
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